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We study the Sinai model for the diffusion of a particle in a one dimensional quenched 
random energy landscape. We consider the particular case of discrete energy landscapes 
made of random ±1 jumps on the semi infinite line Z + with a reflecting wall at the 
origin. We compare the statistical distribution of the successive local minima of the energy 
landscapes, which we derive explicitly, with the dynamical distribution of the position of the 
diffusing particle, which we obtain numerically. At high temperature, the two distributions 
match only in the large time asymptotic regime. At low temperature however, we find even 
at finite times a clear correspondence between the statistical and dynamical distributions, 
with additional interesting oscillatory behaviors. 



09/98 



email: chave@spec.saclay.cea.fr, guitter@spht.saclay.cea.fr 



1. Introduction 



The problem of the diffusion of a particle in a one dimensional quenched random energy 
landscape has been studied for many years ■ If the landscape itself has a random walk 
statistics, many exact predictions exist for the long time behavior of the diffusion process. 
For an unbiased statistics of the potential, i.e in the absence of drift, the average distance 
(x) traveled by the particle grows very slowly at large times, with the Sinai scaling behavior 
(x) ~ (Tln(t)) , where T is the temperature P]. The asymptotic probability distribution 
is universal in the scaling variable x/ (Tln(t)) and its precise form is known exactly for a 
single diffusing particle 0-0]. Moreover, the thermal dispersion (or the distance between 
two diffusing particles in the same random potential) remains finite at large times || . The 
physical picture underlying these results is simply that of a localization of the diffusing 
particle in the deepest energy minimum which it can reach at time t by passing larger and 
larger energy barriers JjH^]. 

These large time asymptotic predictions were tested numerically by computing the 
exact probability distribution for the position x of the particle at finite but large enough 
times t for energy landscapes drawn at random, and then averaging over a large enough 
sample of such landscapes [@-{nj. One should notice that the Sinai scaling (x) ~ (Tln(t)) 2 
defines a very slow diffusive process. For such a process, the long time behavior is expected 
to be reached in practice after extremely long transients, especially at low temperature. 
Such times can be out of reach in practical situations or in finite time numerical simulations. 
This probably explains why the asymptotic predictions are not fully recovered in [Q. In 
any case, very little is known about the finite time behavior of the model or about the 
approach to the asymptotic regime. 

In practice, the localization of the particle in the deepest minimum has the nice 
consequence that some of the dynamical properties of the diffusion process can be deduced 
directly from the corresponding statistical properties of the minima of the random potential. 
More precisely, the underlying idea is to assume that, at times of order t, the particle 
is localized in the deepest minimum of the energy landscape which it could reach by 
passing all the energy barriers AE of height less or equal to V ~ Tln(t). The time 
dependence of various quantities in the dynamical process can then be obtained from the 
equilibrium statistics of the energy minima and its dependence in the highest passable 
energy barrier V. This idea has been implemented recently to develop a new Real Space 
Renormalization Group (RSRG) approach for the Sinai model ||11|| . The RSRG formalism 
consists in performing a suitable decimation of the barriers with height less than V to 
obtain universal Renormalization Group equations for the variation with increasing V of 
the effective energy landscape distribution seen by the particle at the time scale t. Long 
time properties are obtained from the V —>■ oo fixed point of these RSRG equations. This 
technique allows to recover the universal asymptotic probability distribution of a single 
particle for various boundary conditions. It also allows for the prediction of asymptotic 
two-time or two-particle correlations. 

Following the above "statistical picture" leading to the exact asymptotic predictions, 
one may wonder whether this picture could also help to investigate shorter times or very 
low temperatures. In this paper, we test the connection between the dynamical properties 
of the diffusion process and the statistics of the minima of the energy landscapes for finite 
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times. We use a particular distribution of the disorder where the energy landscapes are 
made of successive random ±1 increments, i.e. have the statistics of a discrete random 
walk. The reason for this choice is twofold. On the one hand, we can easily derive explicit 
laws for the statistics of minima at arbitrary (even small) I\ On the other hand, this is 
particularly adapted to numerical simulations of the dynamical process. As we shall see, 
our results corroborate the statistical picture, not only at large V where we recover the 
universal distributions expected for general random energy landscapes with a random walk 
statistics, but also at small V and low temperature, where the discrete nature of the energy 
landscapes is sensible and yields interesting behaviors. 

The paper is organized as follows. In Section 2, we briefly describe the particle diffusion 
process and the distribution for the random energy landscapes. In Section 3, we derive 
explicit formulas for the statistics of energy minima in the presence of energy barriers of 
arbitrary scale I\ Section 4 presents numerical results for the dynamical diffusion process. 
These results are compared with the predictions of the statistical approach of Section 3. 
We gather our conclusions in Section 5. 

2. Diffusion of a particle in a discrete random energy landscape 



Fig.l: A discrete random energy landscape E(x). A reflecting wall prevents 
the particle from exploring the negative x region. The particle jumps to the 
left or to the right with the probabilities given by equation (2.2). 

We consider the following discrete distribution for the quenched energy landscapes. 
On each site x = 0, 1, 2, ... labeled by a non negative integer, we associate a random energy 
variable E(x) such that (see Figure 1) 



The sign ±1 of the increment is drawn at random between successive sites with even 
probability 1/2, i.e. we introduce no drift in the problem. With the choice ( |2.1|) , the 
values of the energy are limited to integers, meaning that we have implicitly fixed an 
underlying energy scale. We can view E(x) as the "height" at site x of a discrete random 
walk describing our energy landscape. We finally set E{— 1) = +oo, i.e. put an infinite 
reflecting wall on the negative side of the origin, while we keep a free boundary condition 



E(x) T 




E(0) = 

AE(x) = E(x + 1) - E(x) = ±1 



(2.1) 
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at x = +00. This choice of boundary condition is clearly not crucial but makes our analysis 
simpler since in this case, deeper minima always appear at increasing values of x. 

For a given realization of the quenched random energy landscape, we consider the 
following discrete time dynamical process: 

- We start at time t = at position x = with energy -E'(O) = 0. 

- Given the position x with energy E(x) at time t, the position at time t + 1 is 

e -f3AE(x) 

x + 1 with probability p x 



x — 1 with probability q x = 



e -pAE(x) , e /3AE(x-l) 

^ (2.2) 



e -(3AE(x) _|_ e f3AE(x-l) 

with AE(x) as in ( |2.1] ) and /3 = 1/T the inverse temperature. 

With the choice ( |2.1| ) for the variation of the energy on neighboring sites, the ratio 
Px/Qx is equal to 1 if the two sites x + 1 and x — 1 have the same energy, or to exp(±2/T) if 
they have different energies. Note that the case x = is special with po = 1 — qo = 1. Note 
also that the particle has to move to one of its neighboring sites at each time step, i.e. we do 
not allow the particle to remain at the same site. As a consequence, the particle occupies 
sites with even position at even times and sites with odd position at odd times. This effect 
results in "residual fluctuations" even at T = 0. The above model is very simple since it 
depends on only one parameter, the temperature T. We study its behavior numerically in 
Section 4 but let us first derive explicit laws for the minima of energy landscapes with the 
discrete random walk statistics (|2.1|). 



3. Statistics of energy minima 

In this section, we forget for a while the dynamical process and focus on the statistics 
of the successive minima of the random energy landscapes which can be reached in the 
presence of increasingly passable energy barriers. Following JT^] , we consider for any given 
realization of the energy landscape the position x = x m i n (r) of the deepest minimum of 
this landscape which can be reached, starting from the wall at the origin x = 0, by passing 
through energy barriers of height less or equal to a fixed scale V. Here x m in and V are non 
negative integers. For each T, we compute the probability pr(x) that x m i n (r) equals x, 
obtained by averaging the quantity 5 X)Xmin (T) over all realizations of the energy landscape. 
Note that, for a given realization, the deepest reachable minimum is not necessarily unique 
but can be degenerate. We will thus consider three different probabilities p^\x), (x) 

and (x), corresponding respectively to the following prescriptions in case of degeneracy: 

(1) we keep only the closest minimum to the origin, 

(2) we keep only the furthest minimum from the origin, 

(3) we take the average over all the degenerate minima, i.e. consider the quantity 
(l/fc)X/i=i^ for degenerate minima at position x^ n . 

' min 

We shall give below explicit formulas for the generating functions 

00 

v¥\z) = J2Pr(*)z X i = l,2,3 (3.1) 

x=0 
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3.1. Random walk in a strip of height V 

Before we proceed, we need to evaluate the probability dr(x) for a random walk 
starting at E(Q) = to have its x first steps inside the strip < E < Y and its (x + l)th 
step at E{x + 1) = —1, i.e. the probability for the walk to leave the strip for the first time 
just after position x and downwards. We recall here the explicit form of the generating 
function Dr(z) = dr(x)z x . In the following, we will always assume that the (x + l)th 
step points downwards, which in practice occurs with probability 1/2, i.e. we will calculate 
2Dy(z) instead of D-p(z). 



r=i 











r-i a 
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Fig. 2: For a strip of height T = 1, the random walk is a series of E = — ► 
1 — > two-step-sequences. To obtain 2Dr(z), each two-step-sequence must 
be dressed with a factor 2Dr-i(z) in the middle. 

The function Dr(z) satisfies the recursion relation: 

Dr(z) ' 



2-z 2 D r -i(z) 



(3.2) 



with the initial condition 2Dq(z) = 1 (for T = 0, the first exit from the zero- width strip 
occurs necessarily just after x = 0). This relation can be understood as follows (see Figure 
2): consider the case T = 1, for which the random walk is necessarily of even length x and 
made of the repetition of x/2 elementary two-step-sequences of the type (E = — > 1 — > 0). 
Each two-step-sequence comes with a factor z 2 in the generating function and occurs with 
a probability (1/2) 2 , leading to, 



°° /r2\2 



xcvcn 



!-(§)' 



(3.3) 



in agreement with ( |3T2"|) and Dq(z) = 1/2. The generating function 2Dr(z) can be obtained 
in the same way by replacing the weight (z/2) 2 for the elementary two-step-sequence by 
a weight (z/2) 2 x 2Dr-i(z), corresponding to inserting between the two steps E = — > 1 
and E = 1 — > an arbitrary sequence of steps limited to the strip of width T — 1,1<E<T 
(see Figure 2). Making this replacement in ( |373| ) directly leads to ( |3.2|) . 
From ( |3.2j ), we can write the generating function Dt(z) as the ratio: 

Pr(z) 



Dr(z) 



Pr+i(z) 



(3.4) 
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where Pr(z) is an even polynomial of degree 2[T/2] in the variable z. Indeed, writing 
D T {z)=P T {z)/Q v {z), (p|) gives 

Pr(z) 1 Qr-i(z) 



Qr(z) 2 -^£=44 2Qr-i(^)-z 2 Pr-iW 

Qr-i(z) 



(3.5) 



leading to recursion relations 

P v {z) = Q V -x{z 



(3.6) 

Qr(z) = 2Q r - 1 (z)-z 2 P r - 1 (z) 



Eliminating Q-p, we get the announced result ( |3.4| ), with, moreover, the recursion relation 

P r +i(z) - 2P T {z) + 2 2 Pr-i(^) = (3.7) 
which determines Pr(z) from the initial values Pq(z) = 1 and P\{z) = 2. Setting 

Pr(z) = z r U r (?) - U r (y) = (f^ Pr (~) (3.8) 

the above recursion relation transforms into 

U r+1 (y) - yU r (y) + U T -i(y) = (3.9) 

with Uo(y) = 1 and U±(y) = y. The solution of this recursion is given by the well known 
Chebishev polynomials Ur(z), characterized by 

sinh((r + l)t) 

U v 2cosh (t) = \\ } } (3.10) 

sinh(ij 

We thus get the explicit form: 



P ( 1 ^ = 1 sinh((r + l)t) 
r Vcosh(t) J (cosh(t)) r sinh(t) 1 1 

and the explicit formula: 

^ ( 1 \ w vSinh((T + l)t) 

D r [z = — — - =cosht — -4t rr 3.12 

V cosh(t) J K sinh((r + 2)t) v ! 

The relation between Dt(z) and the Chebishev polynomials is also derived in |EJ in the 
context of the Temperley-Lieb algebra. In the limit t — > 0, we get the particular result 

Dr(l) = j±± (3-13) 

which is the probability for the walk to leave the strip < E < V downwards at E = 
rather that upwards at E = T. In the limit V —>■ oo, we have 



Dr(z) 1 (3-14) 

which, up to normalizations, is the generating function of Catalan numbers, counting closed 
walks on the semi-infinite line E > 0, as it should. 
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3.2. Explicit forms of 

We now come to the evaluation of the generating functions (z) defined by ( |3.1|) . 
Given an arbitrary landscape of infinite length and for a given V, we have to look for the 
deepest minimum which can be reached by passing barriers of height less or equal to T. To 
do so, we decompose the random walk defining the landscape according to the following 
algorithm: 

1/ Starting at xq = and E(xq) = 0, the walk has a first part entirely included in the 
strip of width r, < E < V. In this part, the last return at E = occurs at some 
position x\ (possibly 0). This defines a first building block of size x\. 

2/ Since x\ is the last return at E = within the strip of size V, two situations may 
occur for the remaining part of the walk: 

2-1/ Either the walk reaches E = T + 1 without returning to E = 0, i.e. the walk 
leaves the strip upwards. In this case, we have reached a barrier too high to be 
passed and we stop the process. 

2-2/ Or E(xi+1) = —1, i.e. the walk leaves the strip downwards. In this case we repeat 
the process 1, i.e starting at position x± + l with E(x± + 1) = —1, we consider the 
following part of the walk entirely in the strip of height T, — l<i?<r — 1 and 
look for the position X2 of the last return to E = — 1 within this strip, defining 
a second block of size x^ — X\ — 1. We repeat the process until it is stopped 
by a too high barrier. The number £ of building blocks is simply related to the 
height 1 — £ of the deepest minimum which could be reached. The position of 
this minimum can be obtained by properly summing the sizes of the blocks. 

This algorithm is illustrated in Figure 3. To compute (z), we simply have to assign 
a weight 2Dr(z) per building block and a factor (z/2) for each descending bond between 
building blocks. The different prescriptions in case of degenerate minima concern only the 
last building block. In the case (1) of the closest minimum to the origin, the size of this 
block does not contribute to the position x m i n of the minimum. This block thus comes with 
a factor 1. Summing over the number £ of blocks, we deduce that V^\z) is proportional 
to (see Figure 3-1): 

P<?> W oc 1 + 2flrWf + (2Z>r(,)|) 2 + . . . = (3.15) 

The normalization factor is fixed by demanding that V^\l) = 1 since, for z = 1, the 
generating function gives the probability to find the minimum at some arbitrary position 
x. We thus get: 

jtf'M = ±^ri\ ( 3 - 16 > 

1 — zDr{z) 

Note that from the previous subsection, we can also understand the prefactor 1 — -Dr(l) as 
the probability for the random walk to eventually leave the strip of size T upwards rather 
than going downwards. For the prescription (2), where we choose the furthest minimum 
from the origin, we have an extra factor 2Dr(z) for the last building block. We thus get 
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Fig. 3: A schematic picture of the decomposition of random walks into build- 
ing blocks. The number £ of building blocks is related to the height 1 — £ of 
the deepest minimum reached without passing a barrier of size greater than 
T. For the different building blocks, we have also indicated the weights ob- 
tained by the average over the random walks to reconstruct the appropriate 

generating function P r . These weights depend on the precise prescription 
for degenerate minima: (1) keep the closest minimum to the origin; (2) keep 
the furthest minimum from the origin; (3) average over the minima. 



in this case (see Figure 3-2): 



V$\z) oc 2D r (z) + 2D v {z)-2D T {z) + (2D r (z) Z -) 2 2D r (z) + ... = (3.17) 
and: 

/r 1 Z) ~ Dr(l) l-zD r (z) 

The difference between the two cases (1) and (2) is thus simply a prefactor Dr(z)/Dr(l) 
coming from the last building block. This factor is nothing but the generating function 
^2 x er(x)z x associated with the probability er(x) to find the furthest minimum at a dis- 
tance x from the closest minimum. 

The prescription (3), where we average over all minima in the case of degeneracy is 
more subtle. Here again, it concerns only the weight of last building block and leads to a 
different prefactor. Suppose we have k + 1 degenerate minima, each minimum will receive 
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a weight l/(k + 1). For the (m + l)'th minimum among the k + 1, we have to assign a 
factor z only to the elementary steps of the walk which are to the left of the minimum 
while the steps which are to the right do not contribute to x^£^ and receive a factor 1 
instead. Still the right part of the block is important to guaranty that the total number 
of minima is precisely k + 1. The generating function between successive minima is thus 
simply 2Dr-i(z)(z/2) 2 on the left side of the chosen minimum and 2Dr-i(l)(l/2) 2 on its 
right side (see Figure 3-3). The prefactor associated with the last block is now proportional 
to: 



i_ j ^ (z*n v _ x {z) 

k+ 1 



^ k + 1 | ^ 

k=o U=o 



1 



z 2 Dr-i(z) _ £> r -i(l) 
2 2 



= 2D r (z)D r (l) 
D r (z)-D r (l] 

where we have used ( [PI) . We finally get 




(3.19) 



V?Hz) = D ^ , In ( x (3.20) 
rK) D r (z)-D r (l) \D r (l)J l-zD r (z) 1 } 

by adjusting the normalization to guaranty that P r (1) = 1. Note that, as for ( |3.18| ), the 
prefactor involves only a function of Dr(z)/Dr(l)- 

The formulas ( |3.16|) - (|3.20|) can be made more explicit by use of ( |3.12| ), leading to 



P <» ( z = _4J\ = x sinh « r+2 w 



cosh(t)/ T + 2 sinh((r + 2)t) -sinh((r + l)t) 



v m ( z = 1 ^ = _L_ x cosh(t) sinh((r + i)t) 



cosh(t) J T + l sinh((r + 2)t) -sinh((r + l)t) 

r + 2cosh(t) sinh((r + l)t) 



p(3) ( z = 1 ^ = cosh(t) sinh((r + l)t) 
r V cosh(t) J sinh((r + 2)t) -sinh((r + l)t) 

sinh((r + 2)t) 

X 



r + l sinh((r + 2)t) 



(T + 2)cosh(t) sinh ((r + l)t) - (T + l)sinh ((T + 2)t) 

(3.21) 

These expressions constitute our central result describing the statistics of minima at finite 
r for the particular distribution of landscapes that we study. 

At large T, an appropriate scaling limit can be obtained by a suitable scaling of z: 

<* <=^+o(p) (3.22) 

This corresponds to setting 

x = T 2 l (3.23) 
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in the original probabilities Pj! (x). In this limit, the three generating functions V r tend 
to the same limiting distribution Voo(s), equal to: 

This distribution is now the Laplace transform 

/>oo 

Poo(s)= / dl Poo (l)e- sl (3.25) 
of the limiting probability density 



o 



1 — >oo 



,(Z) = lmi T 2 pW(x = T 2 l) i = l,2,3 (3.26) 

From (|3.24| ), we get 



oo 

E(2fe + 1)^ 2, 
e-^-H-* ' (3.27) 



The formulas ( |3.24j ) and ( |3.27| ) are actually more general and valid for all the distributions 
of the energy landscape which have a random walk statistics, i.e. tend to a Brownian 
motion distribution in the continuum. These laws can be derived directly by use of the 
RSRG formalism developed in [IT]. They appear there as the large T fixed point of RSRG 



equations describing the flow with T of similar renormalized probability distributions. 
Here, we have made explicit the whole V dependence of these probabilities in the special 
case of an energy landscape made of discrete ±1 increments. 

The nice feature of the formulas Q3.24|) and (|3.27| ) is that they are precisely those 



derived in [|3| for the limiting probability distribution of the dynamical process itself in its 
large time scaling regime. In this case, the natural scaling variable is / = a 2 x/ ln 2 (t), where 
a = 2/T is directly related to the expectation value of ln 2 (p x /l — p x ). We thus recover here 
the asymptotic equivalence between the statistics of minima and the dynamical process, 
with the precise correspondence Y = ln(£)/a. Such a relation, valid in principle only for 
large times, will be tested in Section 4 at shorter times. 

Finally, let us mention that a similar computation of the generating function at finite 
r can be performed for landscapes without wall at the origin. The treatment of the degen- 
eracies is however more involved in this case with even more different possible prescriptions. 
Still, the limiting distribution at large V is independent of the chosen prescription. One 
can also introduce a drift in the problem when drawing the random landscapes. For in- 
stance, choosing in ( |2.1| ) AE(x) = 1 with probability p and AE(x) = —1 with probability 



q = 1 — p, we get for, say, V r ^ the expression 

P« W = '-^'^P, (3.28, 
1 — 2 q z Dy [^4:pqz) 

Different regimes are obtained according to sign of p — q, and to whether p — q is of order 
one (strong bias) or of order 1/r (weak bias). 
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To end this section, let us evaluate the first correction to Voo(s) by further expanding 
the formulas for Pp in powers of 1/r. We get 

_d) ( s\ ( a®\ tanhfv 7 ^) a®( 3 , 2/ /—A , n ( 1 



(3.29) 

with a*- 1 -* = 2, a*- 2 -* = 1 and = 3/2. The above 1/r corrections are probably not univer- 
sal, i.e. they depend on our particular choice for the statistics of landscapes. Still, the 1/r 
correction to the relative differences (V^ — V^)/V^ and (V^ — V^)/V^ should be 

universal. Indeed, they involve the proportionality factors between the different pj, s, i.e. 
the factor D r (z)/D r (l) in ( gig ), or the factor D r (z)/(D r (z) - D T (1)) ln(D r (z) / D r (l)) 
in (|3.20| ). As we already mentioned, D-p{z)/ D-p(l) is the generating function for the prob- 
ability er(x) to have the furthest minimum at a distance x from the closest. The propor- 
tionality factors above thus concern the relative distance between the different degenerate 
minima, and ignore their absolute position. They are the important statistical quantities 
to be used when one is interested in the localization property of the dynamical process. 
If we insist to impose the scaling ( |3.22[ ) appropriate to absolute positions, we get, for 
z = l- s/T 2 



Dr(z) 
D r (l) 



1 + i (l - V2~s coth( v / 2^)) + O (J^j 



D r (z) ^ flMzU =1 If l-v^ooth(y^) \ (1 



D T (z)-D r (i) \D r (i)J r\ 2 / Vr 2 



(3.30) 



At large T, the factors ( |3.30| ) above tend to one. This does not mean that degenerate 
minima disappear in this limit. Indeed, one can easily compute the probability for having 
exactly k + 1 degenerate minima, equal to: 

D r -i(l)\ /£> r -i(l)\ fc T + 2 / T \ k , 

Degeneracies do thus exist even at large T where the above expression tends to (l/2) fc+1 . 
However, most of these degeneracies occur at short distances, and not at distances of 
order T 2 . Relative distances of order T 2 are found only with a probability of order 1/r. 
The 1/r corrections above correspond to situations with exactly two degenerate minima 
at a distance of order T 2 , as corroborated by the fact that averaging over minima gives 
half the value obtained by keeping the furthest minimum. Situations with three or more 
degenerate minima would contribute to higher orders in 1/r. At large T, the scaling 
( p.22|) is therefore not appropriate to deal with the relative position of degenerate minima. 
As noticed in ( |3.14j) , a non trivial large V exists for Dr(z)/Dr(l) without rescaling of z, 
involving the Catalan generating function. Expanding ( |3.14| ) in z, we obtain 

r^oo f ^ x °dd 

e r (x) - e(s) = | _^_^ if x even (3.32) 
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If we admit that relative distances between degenerate minima are the appropriate statis- 
tical quantities to describe the distance between two diffusing particles, this result without 
scaling of x at large times agrees with the idea of localization H] . Still, as explained in |T(J , 
the distribution e(x) behaves like e(x) ~ x~ 3 ^ 2 at large x. Therefore, its moments (x a ) 
diverge for a > 1/2. This divergence comes precisely from the rare configurations (occur- 
ring with probability l/T) with two minima at a distance of order T 2 . They contribute to 
(x a ) by a term (T 2 ) Q /Y = Y 2ot ~ l which diverges at large Y for a > 1/2 (see ||10|| ). 

3.3. Average position of the minimum 

From the formulas for the generating functions (z), it is straightforward to get 

the average position x^(Y) and average squared position x 2 ^ (Y) of the minimum for the 
three prescriptions: 



x«(r) = J2 x Pr\ x ) =z j- z 

x=0 



x 2 (r) = x 2 py i x ) — \ z i~) 

x=0 ^ ( " ' 



vf{z) 

2 = 1 

(3.33) 



Vf{z) 



2 = 1 



Expanding ( |3.21|) in powers of z — 1 around z = 1 (or in powers of t around t = 0), we get 

^ (1) (r) 



2r 2 


+ 5r + 3 




3 


2r 2 


+ 7r + 3 




3 


2r 2 


+ er + 3 



x( 2 \Y) = ^ V (3-34) 
x^(r) 



and 



3 



— =<i) i6r 4 + ssr 3 + i74r 2 + i47r + 45 

x z (Y) = 



,2 



15 

(2). . 16r 4 + 104r 3 + 226r 2 + 179r + 45 
<r- (r) = (3-35) 

—(3), . i6r 4 + 96r 3 + i^r 2 + i63r + 45 

x 2 (Y) = ± 

All of them lead to the same large Y asymptotic formulas 

aj(0(r) r ~°°_r 2 
v ! 3 

( r ) ~ Y5 F 

This asymptotic behavior involves only the scaling variable x/Y 2 and the exact prefactors 
can be directly obtained from the limiting distribution ( |3.24j ) by a suitable expansion in s 
around s = 0. We have in particular the asymptotic relation 

X s ~ — {xf (3.37) 
5 



(3.36) 
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Here again, the difference for the subleading terms in (|3.34| ) and ( 3.35Q for the three 



prescriptions comes from landscapes with two degeneracies separated by a distance of 
order T 2 . Such configurations, occuring with probability ~ 1/r, yield a correction of order 
r 2 /r = r to x(T) and a correction of order T /T = T 3 to x 2 (T). For these situations with 
exactly two effective degenerate minima, it is equivalent to take the average over the minima 
(prescription (3)) and to average over the closest and the furthest minimum (prescriptions 
(1) and (2)). This explains why the subleading coefficient in the prescription (3) is exactly 
the average of the subleading coefficients of prescriptions (1) and (2). Situations with 
three degenerate minima distant from ~ T 2 occur with probability ~ 1/r 2 and influence 
the sub-subleading coefficient. Note also that the constant term in (|3.34| ) and ( |3.35|) must 



be the same for the three prescriptions since, for T = 0, there cannot be any degeneracy 
of the deepest minimum. The above arguments explain why x*- 3 -* (r) is exactly the average 

of x^\T) and x^ 2 \T) and why x 2 (T) differs form the average of x 2 (T) and x 2 (T) 
by a term of order T 2 only. 

In order to have a precise measure of the subleading terms, we introduce the quantity 

(i) = ? 2 "(r) _ ! (3i38) 

(f ) (w">(T)f 

with a factor 12/5 chosen to eliminate the leading large T term, so that Xr ~ 1/r tends 
to zero at large F. Remarkably, in the case (3), an extra cancelation occurs, leading to 
~ 1/r 2 . In the following section, we compare this quantity, as computed from ( |3.34j ) 
and ( |3.35| ), with a similar quantity defined for the dynamical model. 



4. Dynamics 

4-1. Simulation 

Beside the above theoretical statistical predictions, we have made a numerical study 
of the dynamics ( |2.2| ) of a particle for a large sample of quenched random potentials drawn 
to satisfy the relations ( |2 . 1| ) . 

For each drawn landscape, we calculate the probability P(x, t) that the particle sits 
at position x at time t. Due to the discrete nature of the landscape, it is possible to make 
an exact enumeration of all the possible walks arriving at a given position x at some time 
t, and to evaluate their probability deduced from ( |2.2|) for the particular chosen landscape. 

A more efficient way to implement this enumeration is to calculate P(x,t) exactly 
step by step in time from the master equation 

P(x, t + l)= p x -iP(x -l,t) + q x+1 P(x + 1, t) (4.1) 

withp x and q x = l—p x the probabilities of jumping from site x respectively to the right and 
to the left, as explicited in (|2.2|) , with in particular po = 1 — qo = 1 and with the convention 
P-i =0. The master equation is supplemented by the initial condition P(x, 0) = 8 X) q. The 
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computed probability P(x, t) is then averaged over all the energy landscapes of our sample 
(the sample size is 10 5 energy landscapes for most simulations). 

At any finite t, it is clear that the particle cannot reach a position x > t, hence 
P(x,t) = exactly for all x > t. The normalization of the probability requires that 

P(x,t) = 1 at any time t. In practice, the probability to occupy sites far from the 
origin is extremely low. To reduce the computation time, we drop the normalization 
condition and replace it by 

'£ t P{x,t)>l-e (4.2) 

X 

where £< 1. Therefore, instead of describing the complete accessible landscape of size 
x = t at time t, we fix the maximal size at a much smaller value of x such that the 
condition (|4.2|) is fulfilled for all our energy landscapes. Of course, the choice of this size 
depends crucially on the temperature T and on the number of iterations. In practice, we 
take e = 0.01 for most simulations and we checked that this simplification does not lead 
to significant errors. The computation time is significantly reduced for low temperatures, 
since the effectively visited landscape is of much smaller size. With this simplification, we 
were able to study the dynamics up to 10 5 iterations for low temperature regimes (10 for 
T = 2). 

Note that for our particular choice of dynamics, the particle cannot stay at the same 
site for two consecutive times t and t + 1. Thus P(2x + 1, 2t) = P(2x, 2t + 1) = 0. As we 
already mentioned, the net effect of this parity condition is to create residual fluctuations, 
which persist even at T = 0, where barriers of height F = 1 can always be passed. This 
will explain in particular why, at T = 0, some equilibration can take place between all 
the deepest minima accessible by passing V = 1-barriers. In the following subsections, we 
always present results for times of a well defined (even) parity. 

It is instructive to visualize the typical evolution of P(x,t) with time for a fixed 
landscape before making the quenched average over our sample of landscapes. Figure 4 
shows such an evolution with 10 7 iterations. One clearly sees that the regions with a 
large probability of occupation are concentrated around local minima of the potential. 
As t increases, these high density regions migrate to deeper minima. The duration of 
occupation of a local minimum in logarithmic scale is roughly proportional to the height 
of the energy barrier on its right. 



4-2. Distribution P(x,t) 

Let us first present our numerical results for the distribution P(x,t). We will use 
overbars to denote the average over our sample of landscapes, while brackets (■) will denote 
the thermal average estimated from the probability P(x, t) computed for a fixed landscape. 
We are interested here in the average distribution P(x, t) at a, fixedlaxge enough time. From 



the asymptotic large V results (|3.26|) and (|3.36|) , we expect that, at large t 



(x)(t)P(x,t) 




(4.3) 



with poo given by ( |3.27| ) . 
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Fig. 4: Evolution with time (in logarithmic scale) of the distribution P(x,t) 
in a given energy landscape (drawn below). The evolution runs over 10 7 
iterations. The intensity in the greyscale is proportional to — In P(x,t), i.e. 
darker regions correspond to higher values of P(x,t). 

Figure 5 shows our results for T = 1/2 and different values of t. The agreement 
with the asymptotic exact formula is apparently very good. To have a better quantitative 
evaluation of how close we are to the asymptotic result, we will study in the next section 
the first and second moments of the distribution. As we shall see, significant deviations 
do actually exist, some of which can be well explained by our finite V corrections to 
the asymptotic statistics. We will also discover some interesting underlying oscillatory 
behaviors. 

4-3. Results for the first and second moments and comparison with the statistics of minima 

We present here our numerical results for the average position (x)(t) and average 
squared position (x 2 )(t) for varying time t and at various temperatures. We first check if 
(x 2 )(t) and ((x)(t)) 2 obey the expected asymptotic relation ( |3.37| ) with a proportionality 
factor 12/5. Note that for a purely diffusive system in a homogeneous medium (i.e. in 
a flat energy landscape) where (x)(t) ~ \ft 1 a similar relation holds, but with a smaller 
proportionality factor n/2. 
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Fig. 5: The rescaled average distribution (x)(t)P(x,t) as a function of the 

rescaled variable x/(x)(t) at T = 1/2 and for various values of the time 
t = 1 000, 5 000, 10 000, 50 000 and 100 000 The solid line indicates the exact 
asymptotic formula as given by ([Of) and (|3.27[). 



Figure 6 presents the corresponding data for several temperatures from T = 2/5 to 

T = 2. For the range of (x) presented here, the asymptotic formula ( |3.37| ) and the 
more complete statistical relations obtained from our finite V predictions by eliminating V 
between ( |3.34| ) and ( |3.35| ) do not differ significantly. We therefore expect, if the dynamics 
follows the statistics of minima, that the asymptotic linear relation is verified in the whole 

2 

range of (x) . This is precisely what we observe at low temperatures (below T ~ 1). In this 
regime, a more refined analysis will reveal very interesting underlying oscillatory behaviors, 
as emphasized below. At high temperatures (above T ~ 1), we see a significant deviation 

2 

from the expected law, with a behavior closer to a purely diffusive regime at small (x) . 
Still, even at T = 2 (see inset in Figure 6), the correct slope is eventually recovered at 

2 

large (x) . In this range of temperatures, the observed deviation is not explained by our 
finite r corrections to the asymptotic law p^, but is more simply the effect of a short 
time diffusion regime in which the particle does not yet feel the random potential. Such a 
regime will last until the particle reaches a distance x such that E(x) ~ y/x ~ T, i.e up to 
a transition time t ~ x 2 ~ T 4 . 

For a comparison of our data with the statistics of minima of the previous section, we 
need to consider instead of the first and second moments more refined quantities which in 
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Fig. 6: A numerical check of the relation ( |3.37| ). The solid line has the slope 
12/5 expected from ( |3.37|) while the dashed line has the slope 7r/2 expected 
for the normal diffusion is a flat landscape. For T = 2, a broader window of 

2 

(x) (see inset) is needed to recover the asymptotic law. 

practice contain exactly the same information but are more adequate for our purposes since 
they clearly emphasize the finite F corrections. Anticipating our conclusions, we focus on 
the statistics (3) which corresponds to an averaging over degenerate minima. Inverting 
the formula ( |3.34j) for the prescription (3), we consider instead of (x)(t) the equivalent 
quantity: 



-3 + \ 3 + 6(x)(t) 
F(t) = V 2 (4.4) 

such that x~( 3 \r(t)) = (x)(t). The quantity T(t) is thus an estimate of the effective height of 
the barriers which can be passed at the time scale t, obtained by matching the first moment 
measured in the dynamical process with the average position of the minima resulting from 
the formula (|3.34| ) in case (3). 



Figure 7 shows T(t) as a function of ln(t) for different values of the temperature. At 
low temperature, F(t) oscillates around an average straight line and develops plateaus at 
integer values of V. These plateaus are of course a signature of the underlying discrete 
nature of the landscape, and are an indication of the actual relation between the dynamics 
and a process of passing increasing discrete barriers. This effect disappears at higher 
temperatures. Discarding these oscillations, the curves have an average linear growth with 
a slope directly proportional to the temperature. More precisely, we can reasonably fit 
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Fig.7: The quantity T(t) as defined by Eq. (|Oj) for T = 1/20 to T = 2, as 
a function of ln(£). The solid lines show the corresponding fits of Eq. (4.4) 
for the average linear growth of the curves remaining after discarding the 
superimposed oscillations. The use of the complete formula ( |4.4j ) is necessary 
to obtain the correct position of the fit at low temperatures. At high tem- 
perature (T = 2 in the inset), the fit (4.4) is better if we define T(t) by the 



asymptotic formula T(t) ~ y (3/2){x)(t) 

this average linear growth by the formula (see Figure 7: for the oscillating curves, the fits 
reasonably match the maxima of the oscillations) 

T(t) = ^ln (f) + 1 (4.5) 



2 V*o 

with ln(io) ~ 4. The proportionality factor 1/2 = T/a between T(t) and Tln(t) is that 
expected from the correspondence between statistics and dynamics in the asymptotic limit, 
as already discussed. At low temperature and shorter times, this correspondence is still 
reasonably good, apart from the superimposed oscillations. The additive constant 1 in 
( |4.5|) can be understood as the effect of the "residual fluctuations" which remain at T = 
from the parity condition and make the barriers of height T = 1 always passable. To obtain 
the correct position of the linear fit, we definitely had to use the complete formula 



which presents a shift of —3/2 with respect to the asymptotic relation T(t) ~ y (3/2)(x)(t) 



obtained from ( 3.36 ). This reflects the importance of the finite V corrections at low tem- 
peratures. For higher temperature however (T = 2), we obtain a better scaling with the 
asymptotic law without shift (see the inset in Figure 7, black triangles) than with the shift 
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(white triangles). Again at high temperatures, the deviation from the asymptotic limit is 
not explained by finite V corrections alone. 

We now analyze our data for the second moment (x 2 )(t). Here again, we prefer to 
consider the more adequate quantity 



x(t) = _ , (46) 



(¥) (MM) 

copied from ( |3.38j ) to measure the deviation from the asymptotic regime ( |3.37| ). 




Fig. 8: The quantity x(t) defined by Eq. (\i.6\j as a function of T(t) for tem- 
peratures ranging from T = 1/5 to T = 2. The solid line shows the statistical 

estimate x[?^(r) as defined by Eq. ( |3.38|) . The cross in the inset indicates 
the small corrections at T = 1 to this statistical value due to parity effects, 
as computed in Appendix A. The dashed line in the lower inset indicates the 
value of x f° r pure diffusion in flat landscape. 

Figure 8 presents our results for x{t) versus and a comparison with the corre- 

(3) 

sponding statistical relation calculated from the previous section between Xr an d T. We 
recover again a low temperature regime with oscillations and a high temperature regime 
without oscillations. In this high temperature regime, the asymptotic limit is reached very 
slowly and the short time dependence is different from the finite V predictions for the 
statistics of minima. The subleading short time corrections to the universal asymptotic 
behavior are thus different from that predicted by the statistics of minima. At low tem- 
perature however, the oscillations of x are peaked around integer values of T, and sharpen 
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as the temperature decreases to zero. In this limit T — > 0, the height of the peaks tend 

(3) 

precisely to the value of Xr a ^ the corresponding integer. We thus recover our predictions 
for the statistics of minima, which strictly speaking are valid only for integer values of the 
barrier V. We also observe two additional dynamical effects: 

- Between two consecutive integer values V — 1 and T, we find a transition regime 
with a strong depletion of x- We interpret this effect as resulting from a period of 
equilibration of the particle passing from the statistics of minima at scale T — 1 to 
that at scale V. 

- As the temperature increases, the peaks are rounded and their maxima slowly move 
to lower values of x- The peaks eventually disappear at high temperature. 

We interpret the above results as follows. From the formula (|4.5| ), the barriers of 
height r are passed at times of order t{Y) ~ i ex P[2(r — 1)/T]. After passing these 
barriers, we admit that the time needed for equilibration in the (always present) degenerate 
minima is itself of the order of a multiple of t(T). As the temperature is lowered, the 
corresponding proportionality factor remains finite due to the residual fluctuations. During 
the equilibration process at a given T, the data get closer to the equilibrium distribution 
of the minima. The time needed to pass the next barriers, i.e. those of height T + 1 is 
from ( |4.5| ) t(T + 1) = t(T) exp(2/T), that is again a finite multiple of t(T), depending only 
on the temperature. In the low temperature regime, the particle thus has enough time to 
equilibrate and recover the statistics of minima for a fixed passable height V before it starts 
passing the barriers of height T + 1. Conversely, at high temperature, the particle keeps 
finding better and better minima by passing increasing barriers without equilibration for 
each passed V. In particular, the particle does not feel the discrete nature of the potential. 
This results in the suppression of the oscillations and a behavior closer to normal diffusion. 

In the above analysis, the use of statistics (3) is crucial on the one hand to get peaks 
precisely at integer values of V(t) and on the other hand to recover the theoretical value of 
x(r) at the peak for low temperatures. These conditions eliminate the two other statistics. 
At this level of precision, the purely asymptotic result x = is also ruled out. As discussed 
in the Appendix A, the statistics (3) itself must be modified by very small corrections due 
to parity effects, i.e. the fact that the particle at time t cannot sit right at the correct 
minima if those happen to have a parity different from t. This effect is sensible only for the 
peak at r = 1. The corresponding correction is calculated in the Appendix A and leads to 
a very small reduction (~ 3.5%) of the peak (see Figure 8), consistent with the data of the 
dynamical process. The parity correction is more important for quantities which measure 
the localization of the particle, as we discuss in the next section. The use of a different 
dynamics allowing the particle to remain at the same site should suppress this correction. 

4-4- Localization 

In order to have an idea of how localized is the particle, we have measured the prob- 
ability second moment: 



which estimates the probability that two independent particles evolving in the same 
quenched potential arrive exactly at the same site at time t. This quantity is similar 




(4.7) 



X 
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to the participation ratio for localized quantum particles [[13]]. A non zero value of Y 2 is 
the signal of a localization. We must however distinguish between two different effects 
which run counter to this localization and therefore lower Y 2 : the effect of temperature 
which broadens the distribution of a particle around its average position in a minimum, 
and the existence of several degenerate minima in which the particles can fall. It is this 
last effect that we measure at very low temperatures. If we assume that, when T — > 0, the 
particle is localized exactly in the deepest minima at scale T(t), we estimate Y2(t) by 



Y (D = V 1 F + 2 ( F ^ ' 

21 ] A^ + i2(r + i) Ur + i)J 



r + 2 n /2(r + i) 

in 



r + 2 



with r = T(t). In the formula above, we used the probability ( |3.31| ) to have exactly k + 1 
degenerate minima, weighed by the probability l/(/c + 1) to have the two particles in the 



same minimum. In practice, as we already noticed, the formula (4^) must be corrected to 
account for the parity effects of the dynamics. At large T, this simply amounts to reduce 
Y 2 by afactor (r + 2)/(2r + 3) + (l/2) x (T + l)/(2r + 3) since Y 2 must typically be divided 
by 2 if the minima have the wrong parity (which occurs with probability (T + 1)/ (21? + 3) 
for even times, see Appendix B) and the particle has to sit on the two neighbors of the 
minimum. With this estimate, Y 2 (t) thus tends at large times (T(t) — > oo) to a finite 
value (1/2) ln2. The same formula ( |4.8|) with the above reduction factor gives also a good 
estimate at small T, with for instance a value Y% ~ (12/5) ln(4/3) ~ 0.69 at T = 1. In 
Appendix B, we present a detailed analysis of the exact corrections for the statistical model 
at T = 1, leading to a value Y 2 ~ 0.706415. 

Figure 9 shows our numerical results for Y 2 [t) as a function of T(t). At low tem- 
peratures, we recover peaks at integer values of V. At low temperatures, the value of 
the peak at V = 1 is in perfect agreement with the estimate of Appendix B and our re- 
sults are consistent with a localization of the particle in all the degenerate minima. At 
high temperatures, Y 2 (t) decreases with time but still tends to a finite value at large T, 
apparently proportional to 1/T. We have also measured the more usual Renyi entropy 
H 2 {t) = — ln(^ x P 2 (x, t)) 0. In contrast with Y 2 , which is an average value over the 
disorder, exp[— H 2 ] gives the typical value of the probability for two particles to be at the 
same site at time t. We find that exp[— H 2 ] and Y 2 display the same behavior and differ 
by a roughly constant multiplicative factor. 

Another quantity of interest for the measure of the localization of a particle is the 
dispersion, defined by 

Ax 2 (t) = (x 2 )(t)-((x)(t)) 2 (4.9) 

If the particle were localized in a single minimum, the dispersion Ax 2 (t) would not grow in- 
definitely with time (or with T) but rather would reach a finite limit. The situation is quite 
different if, as we expect, the particle is localized in several degenerate minima since, as we 
already mentioned, minima can be separated by a distance of order T 2 with a probability 
1/r. In this case, Ax 2 grows like (r 2 ) 2 /r = T 3 and is thus infinite asymptotically. 
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Fig. 9: The quantity Y 2 {i) as denned in Eq. (|4.7| ) as a function of T(t) for 
temperatures ranging from T = 1/5 to T = 2. The horizontal bar is the 
value of the peak estimated in Appendix B from the statistics with parity 
corrections. 



Figure 10 shows our numerical results for Ax 2 (t) as a function of F(t) in logarithmic 
scales. At low temperatures, the peaks of Ax 2 (t) follow the statistics (3) of the degenerate 
minima, with the additional parity corrections, as computed for V = 1 in the appendix B. 
We thus recover a regime of localization in all the degenerate minima. At high tempera- 
tures, Ax 2 (t) increases rapidly with time, with a scaling compatible with the expected T 3 
dependence. 



5. Conclusions and discussion 

We have studied the problem of the one dimensional diffusion of a particle in a semi- 
infinite quenched random energy landscape with discrete integer heights and for different 
temperatures T. Our data converge to the expected large time universal asymptotic limit, 
but very slowly (as l/ln£) and the approach to this limit crucially depends on temperature. 
To quantify the finite time corrections, we have compared our numerical data for the 
diffusion process at time t with exact results for the statistics of the local minima of the 
landscapes reached by passing energy barriers of increasing (integer) size T. Below T ~ 1, 
we find a low temperature regime in which the dynamical process precisely follows the 
statistics of minima which corresponds to average over all the degenerate minima at a 
given scale V with the expected correspondence r <-> Tint. Our data are consistent with a 
localization of the particle equally distributed in all these degenerate minima. Interesting 
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Fig. 10: The quantity Ax 2 (t) (in logarithmic scale) as defined in Eq. (|4.9| ) 
as a function of F(t) (in logarithmic scale) for temperatures ranging from 
T = 1/10 to T = 2. The horizontal bar indicates the value of the peak 
estimated in Appendix B from the statistics with parity corrections. 

transition regimes interpolate between integer values of the scale V. At higher temperatures 
above T ~ 1, the approach to the asymptotic laws do not follow the finite Y corrections to 
the statistical laws, at least in the regime of times that we consider in our simulation. 

We interpret these two different regimes as follows: at small T, the particle has enough 
time to equilibrate and recovers the statistics at a given Y before it starts passing barriers 
of heights Y + 1. At high temperatures however, the particle keeps passing barriers of 
increasing Y without having enough time for equilibration at a fixed Y. 

A natural question is whether the existence of two regimes is only a finite time effect, 
or whether it persists for larger times. In the first case, the correspondence between 
the dynamics and the statistics should be recovered at larger times. In the second case, 
it is tempting to expect a transition temperature between the two regimes. The low 
temperature regime would be anyway a consequence of an underlying discrete cut-off for 
the steps of the energy landscape. The transition temperature would thus be of the order of 
this cut-off and would tend to zero in the continuum limit where only the high temperature 
regime would persist. 
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Appendix A. Corrections due to residual fluctuations at T = 

We will concentrate here on our statistics (3), which consists in averaging over all the 
degenerate minima with equal probability, that is, in the case of k degenerate minima, in 
assigning a weight factor 1/k to each minimum. This statistics implicitly assumes that, 
after a transition period of equilibration, each degenerate minimum is visited with even 
probability. 

For our particular choice of dynamics however, the particle is not allowed to remain at 
the same site for two consecutive times t and t+1. As we already mentioned, this results 
in a parity condition, namely that the particle occupies even sites at even times and odd 
sites at odd times. Thus, even at T = 0, the particle has to fluctuate from odd to even 
sites. This is what we called "residual fluctuations" , which in particular make barriers of 
height T = 1 always passable, even at T = 0, eventually leading to some equilibration 
between all the deepest minima accessible by passing these V = 1-barriers. 

Another effect of these residual fluctuations is that, by preventing the particle to re- 
main seated right in the minima for all times, they create small corrections to the estimates 
of statistics (3), in particular for the limiting values of the peaks at T — > of the quantities 
X, Y 2 and Ax 2 . 

As we will now discuss, these corrections are of two types: a parity correction and a 
wall correction. 

Parity correction: Let us consider a time t with a given parity, say even, so that the 
particle is forced to occupy sites with the same parity. For a fixed landscape and a fixed 
r, it is clear that all the degenerate deepest minima have a well defined, common parity 
since returning to the same height requires an even number of steps. This parity however 
may or may not be that of t. If the two parities coincide, then the particle can fall at 
T — > precisely in the minima and our calculations using statistics (3) are valid. If 
the two parities do not match, the particle cannot fall right in the deepest minima, but 
will rather occupy the two neighboring positions on both sides of each minimum. In 
this case, rather than averaging the positions x^^, m = 1, . . . , k of the k minima with 

weights 1/k, we should average the positions x^l ± 1 with weights l/(2k) (of course, if 
two minima are distant by two elementary steps only, the point in between should receive 
a weight l/(2/c) + l/(2fc) = 1/k). What is the probability for the minima to have the 
wrong parity? Looking at, say, the closest minimum, we can select the even minima by 
considering the combination (V^\z) J r T > ^\—z))/2 and the odd minima by considering 
(V^ \z) — Vy \—z))/2 instead. The probability for a minimum to be even or odd is thus 



As far as x is concerned, it is clear the parity shift creates no correction since x ™l is 
precisely the average of x^l — 1 and x^l + 1. For x 2 however, it results in a shift by a 




(A.l) 
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factor: 

E Yk fa™ - 1)2 + ^ + 1)2 ) - l {x ^ 2 = 1 ( A - 2) 

m=l 

irrespectively of the number of degenerate minima. For even times, we thus get a parity 
correction to x 2 equal to 1 x (T + l)/(2T + 3). This correction is negligible at large V but 

can be measured at T = 1, where it predicts a shift of 2/5 to a; 2 (1). The same effect is 
more sensible for Y<i and Ax 2 , and will be discussed in the Appendix B. 
Wall correction: Another correction comes from the presence of the wall at x = in the 
particular case where x m i n = itself is a minimum. This situation occurs with probability: 

^(0) = ^ (A.3) 

Let us thus assume that x^ in = is the first deepest minimum, together with (k — 1) other 

minima at even positions xj^l, m = 2, . . .k. Let us also assume that t is even so that 
the particle can sit precisely in these minima. Still, after equilibration, the minimum at 
^min = is less probable that the other minima. This effect is visible at T = and T = 1 
where the equilibration is due to residual fluctuations only. 



Weights at odd times 




Weights at even times 



Fig. 11: The modified equilibrium weights for k degenerate minima at V = 1 
in the case of a first minimum right at the wall. Each weight at even times is 
the average of the neighboring weights at odd times and conversely. 

In this case (see Figure 11), the k minima are at positions x^ n = 2m — 2. Since the 

minimum at x^ n has no accessible neighbor on its left, it is easy to see that, at even times, 
this minimum is less probable that the others by a factor of 2, leading to a probability 
l/(2k — 1) for this minimum and a probability 2/ (2k — 1) for the (k — 1) others, instead of 
an equal probability 1/k for each minimum. The correction to x^(l) is thus: 

( 2 1\ A . k- 1 . A , 

v 7 m=2 
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(3) 

and that to x 2 (1) 



(s^T-s) = ft*- 1 ) < A - 5 > 

v 7 m=2 



According to (|3.31|) and (| A.3| ), such a situation occurs with probability (1/3) x 
(3/4X1/4)*- 1 

Combining the parity correction and the wall correction, we get 

(3)r^ 1 >^3 /l^" 1 92-31n(3) , . 

fc=i v 7 

instead of x^ 3 \l) = 11/3, i.e. numerically 3.696 instead of 3.666, and 

-5 -5(3), N 2 1 ^3/l\ fc " 1 2 /I . 4724 . . . 

k=l V 7 



(3) 

instead of x 2 (1) = 932/27, i.e. numerically 34.99 instead of 34.52. With these values, 
the estimate = 0.06978 is modified into x C orr. = 0.06733, i.e. is lowered by 3.5% only. 



Appendix B. Computation of Y2 and Ax 2 in the limit T — > 

In this Appendix, we will calculate the value of the peak at T — > and r = 1 of Y2 
and Ax 2 . These value are very sensitive to the residual fluctuations. We will thus follow 
the same lines as in the previous Appendix. We will consider even times and a situation 
with k degenerate minima at positions x^ n = x^ n + 2m — 2, m = 1, . . . , k. According to 
( P-31|) , this situation occurs with probability (3/4)(l/4) fc_1 . We distinguish between three 
possibilities for the first minimum: 

(i) x^J n = 0. According to ( | A.3|) , this occurs with probability 1/3. 

(ii) cc£j n is odd. According to ( | A.l[ ),this occurs with probability 2/5. 

(iii) x£j n is even and non zero. This occurs with probability 1 — 1/3 — 2/5 = 4/15. 

The weights associated to the minima (or their neighbors) in the three situations are 
depicted in Figure 12. The corresponding contributions to Y2 are 

1 \ 2 „ 2 \ 2 4A;-3 

+ ( k ~ !) 771 T = 1T2 in case W 



2k-\) v ' \2k-\) (2k- l) 2 
1 \ 2 „ fl\ 2 ( 1 \ 2 2k -1 



21/ +( *- 1 HfcJ + UJ ="2P- incase(ii) (B ' i; 

n 2 1 



k ( — J = — in case (iii) 
k I k 
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Fig. 12: Respective weights of accessible minima at even times for V = 1 in 
the three situations: (i) the first minimum is right at the wall; (ii) the first 
minimum is at an odd position; (iii) the first minimum is at an even position 
but not at the wall. 

Combining these contributions with their respective weights, we get 

fc-i 



1 4k -3 2 2k - 1 4 1 
4 \4j \3 X (2k- l) 2 + 5 X 2k 2 + 15 X kj ( B .2) 



k=i 
0.706415 



Similarly, the contributions to Ax 2 are 



2 



v-/ , 2 / 2 \ 4k(k 3 - 2k 2 + 2k - 1) 

(B.3) 



k 

^(*£L " + 1 E + 2 ™ " 3 ) 2 + ^(^„ + 2fc - I) 2 

m=2 



f^(^n-l) + ^E(^n + 2m-3) + ^(4 1 L + 2/ C -l)') = in case (ii) 



(B.4) 
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k / k \ 2 

£ X>£L + 2 ™ - 2 ) 2 - ^E(^in + 2m-2)j =*-^i incase(iii) ( B.5) 

m=l \ m=l / 

Combining these contributions with their respective weights, we get 

a2 ^3/l\ fc_1 fl 4k(k 3 - 2k 2 + 2k - 1) 2 k 2 +2 4 A; 2 - 1 1 

Ax 2 = > - - < - x — — + - x 1 x } 

^4V4y \3 3(2*; -l) 2 5 3 15 3 J 

= 0.789288 

(B.6) 
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